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Abstract 

We examine the quark mass dependence of the pion vector form factor, particularly the curvature (mean quartic 
radius). We focus our study on the consequences of assuming that the coupling constant of the p to pions, gp-mr, 
is largely independent of the quark mass while the quark mass dependence of the p-mass is given by recent lattice 
data. By employing the Omnes representation we can provide a very clean estimate for a certain combination of 
the curvature and the square radius, whose quark mass dependence could be determined from lattice computations. 
This study provides an independent access to the quark mass dependence of the pnn coupling and in this way a 
non-trivial check of the systematics of chiral extrapolations. We also provide an improved value for the curvature for 
physical values for the quark masses, namely (r") = 0.73 ± 0.09 fm"* or equivalently cv = 4.00 ± 0.50 GeV"*. 
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1. Introduction and notation 



The expectation value of the vector current between two pion fields may be written as 

{n^{p')\r\7r^{p))^{p + prFnq') ■ 

Since the only form factor that we discuss is the charged pion form factor, we will denote it simply as F{q^), 
where q — p ~ p' . It in conventionally normalized as F[Q) = 1. 

An expansion around zero momentum transfer allows for a physical interpretation of the form factor in 
terms of the pion's rest frame charge density p{r), given by 



F{t) = l + -{r^),t+-{r^),t' + 0{t') . (1) 



Here we used the notation of Ref. [1 . Alternatively, in Ref. [5] the curvature of the form factor was introduced 
via 

F{t) = l + l{r^)t + c^yt' + 0{t^) . (2) 

D 



Comparing with Eq. ([T]) we see that 



The first term in the form-factor expansion is the conventional charge normalization J (Prp{r) = (!) p = 1, 
and the derivative at the origin provides the (vector or charge) pion radius (r^) = &{dF / dq^){Q). At the next 
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Table 1 

Current theoretical estimates of the pion's quartic radius . The lattice results are our estimate based on the form factor data 
from |15l . 





Method 


3.3 


VMD [2] 


2± 2.5 


Lattice m,r = 0.33 GeV [E] 


2± 2.3 


Lattice extrapolated to m,r = 0.139 GeV |15| 


4± 2 


NNLO xPT [16] 


3.1 ± 0.4 


Padc approximants |17) 


3.6 ±0.6 


Eq. (0 this work 



order, the non-relativistic interpretation should be modified as effects of boosting the pion wave function 
should begin to appear. In this article we will ignore this subtlety, and simply use equivalently the term 
pion's mean quartic radius or form factor curvature. We will mainly focus on this quantity. Using both chiral 
perturbation theory and dispersion relations we find a reliable value for the curvature. 

Assuming that the ptttt coupling is largely independent of the quark masses for a given quark mass depen- 
dence of TOp, we can also predict the quark mass dependence of the curvature. The quark mass dependence 
of the p properties was studied in various recent lattice simulations [3l4j as well as using unitarized chiral 
perturbation theory [5]. 

Abundant data on the pion form factor exists, that can be obtained from the Durham reaction database. 
For the timelike form factor we use contemporary sets from the CMD2, KLOE, and SND experiments 
[6l7l8j . In addition there is higher energy data from Babar[^ that shows the p(1700) and a shoulder that 
could correspond to the /5(1400). However, in our study we employ only the p(770). We will therefore not 
extend our study beyond 1.2 GeV, where in addition KK and other inelastic channels start to contribute 
significantly. It is this condition that prevents us from studying the radius instead of the curvature, as will 
be explained below. 

In the case of the spacelike pion form factor the data is taken from the CERN NA7 collaboration [10 . 
The more recent data from JLAB [TT] was taken at values of too large for our study. For a recent review 
on the status of the spacelike form factor see Ref. . 

The two-loop chiral perturbation theory (^PT) analysis of [13] yielded a mean quadratic radius of 



(r^) = 0.452(13) fm^ 



(3) 



which is the currently accepted value [14j. In order to get a feeling on what to expect for the quartic radius, 
we start with some simple classical examples. For this discussion we will divide it by the mean square radius 
squared, the resulting ratio R = (r^)/(r^)^ quantifies the radial spread of the charge distribution. For a 
charge conducting sphere the spread is minimal with R = 1 (all the charge is at the surface), and for a 
uniformly charged dielectric sphere R — 25/21. On the other hand, the ratio is as large as 5/2 for a charge 
distribution with an exponential dependence on the radius e^™''. A vector-meson pole form factor 

gives an even higher value of 10/3 [2]. 

Some results about the pion's quartic radius, including this work, is collected in table [TJ Furthermore, 
after having analyzed the pion em form factor data by using analyticity, the curvature was constrained in the 
range [0.25 GeV"^,7.57 GeV"''] in Ref. 18J and [2.3 GeV~^,5.4 GeV"^] in a very recent analysis [19|. Of 
particular interest for us, and for a lattice determination, is the quark mass dependence, or dependence, 
of the quartic radius. A study of this within chiral perturbation theory would require control of the N'^LO 
Lagrangian, since the quartic radius is NNLO itself, and this seems out of today's reach. 

We examine the problem with the simphfying assumption that .^pt^tt 

is m^-independent while the 
dependence of nip is taken from other sources. To control the model dependence, we employ the Omnes 
representation of the form factor, sketched in Subsection 12.11 below . In the absence of form factor zeroes, and 



2 



neglecting inelastic channels, this only requires knowledge of the elastic pion-pion scattering phase shifts. 
We parametrize them, with a simple Breit-Wigner model described in Subsection 12.31 Since this model 
contains gpTnr and nip as the only parameters, the above assumption can be employed in a straightforward 
way. In addition we also use the predictions of unitarized chiral perturbation theory for both rho mass and 
coupling. The resulting quark mass dependence of the rho properties were investigated in Ref. [5]. Using 
this alternative parametrization we get almost identical results. The pion mass dependence of the curvature 
turns out to be similar to that of the square radius. 



2. Omnes representation of the form factor 



2.1. Basics 



The Omnes equation encodes the analyticity properties of the pion form factor F{s), that has an elastic 
unitarity cut on the positive s-axis for s G (4to^,oo), and is otherwise analytic. Further superimposed cuts 
due to inelastic channels are neglected in its derivation, and the form factor is assumed to have no zeroes 
(which, as we know today, is phenomcnologically correct). Wc have explored the possibility of zeroes in 
the complex plane by analytically continuing the experimental data with the help of the Cauchy-Riemann 
equations [H]. For a small band around the real axis, they can be excluded. Some remarks on inelastic 
channels can be found in 

The starting point is the well known relation Im(i^) = tan(5iiRe(-F), which relates the discontinuity in 
the vector form factor to the elastic scattering phase shift in the vector-isovector channel. From this relation 
the Watson theorem follows straightforwardly. Since the large-g^ asymptotic behavior of the form factor is 
known from QCD counting rules [23,, F{c[^) a matter of principle one may write an unsubtracted 

dispersion relation, which reads for arbitrary t 

00 

F{t)^- f dstanSnis)^^^^^ . (4) 
n J s — t — le 

We specified "as a matter of principle" since the QCD counting rules apply only when elastic scattering is 
irrelevant by the numerous inelastic channels open. However, in this work we only want to use low energy 
input (up to 1.2 GeV) and we will therefore use a subtracted dispersion relation below and cut the high 
energy contributions with a cut-off. The variation of the results with this cut-off provides a systematic 
uncertainty in our work, which, as a consequence of the subtraction, turns out to be moderate. 

If there are no bound state poles, as is the case of tttt scattering for physical quark masses, nor the form 
factor vanishes anywhere in the complex plane, as we presume for F{t), the celebrated solution family of 
this equation provides a representation of the form-factor in terms of the scattering phase, known as the 
Omnes representation. The standard treatment proceeds by deriving an integral equation for \ogF{t)/{2i) 
instead of F{t) itself, 

4m2 4m^ 

Instead of this relation we may use a subtracted version. This will allow us to effectively suppress the high 
energy behaviour of the phase shifts. In particular we will use a twice subtracted version which reads after 
exponentiation 



F(t) = exp 



/ 00 \ 

P.t+'- fds 



TT J s'^{s — t — ie) 

4)nJ J 



(6) 
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Note, the normalization condition of the form factor prohibits a constant term in the exponent. The constant 
Pi can be identified with the square radius of the pion 

Pi = {r^)/6 . 

This representation of the form factor has been used in the hterature, see for example [24j . 

2.2. Phase-shift representation of the quartic radius 

Recalling the definition of the curvature of the pion form factor (c.f. Eq. ([2])) we may read off an expression 
for Cy directly from Eq. ([6|): 

oo 



120 72^^ ' ■ TT 



ds"-^ (7) 



4m; 



which is quite a beautiful formula, since it allows a third independent extraction of the curvature Cy besides 
NNLO xPT or a fit to spacclike data beyond the linear fall in t where uncertainties get large. Instead we 
employ only the elastic phase shift. In addition, since the quantity 

cl^cl-j^{rY (8) 

is described solely in terms of the tttt p-wave phase shifts, its quark mass dependence is closely linked to that 
of the p-meson properties. This relation is analogous to others existing for the mean-square radius [2l25j . 

2.3. A simple Breit-Wigner model 

To provide an estimate of the form factor based on the Omnes representation, we employ a simple rel- 
ativistic Breit-Wigner model of the scattering amplitude, in which an s-channel resonance dominates the 
scattering 

c 

"-iii^) = 2 ■ r TT 

s-mj- impVtotys) 

where Ftot is the total width of the p resonance with 



,2 „3 „2 ^2)3/2 



y pTTTl 



(10) 



We neglected terms of order Fj^j . Here c is an irrelevant constant that may be expressed in terms of nip and 
Etot- We may write 

Sii{s) = arctan^^— — arctan _ fii\ 

Reaii(s) nip — s 

With this phase variation the integral representation converges without additional subtractions (although 
the high-energy tail is ad-hoc), however, even values as large as s = 7 GeV^ contribute to the integral. 

A good fit can be seen in Fig. [T] for the phase shift and square form-factor modulus. To produce the 
figures we use mf'y^ = 139 MeV, nip = 775 MeV, Tf'y'' = 150 MeV (to determine gp^^). 

Using the formula given above we may now extract the curvature directly from the elastic pion phase shifts, 
using the square radius as input. We find 

Cy = 3.75 ± 0.33 GeV^'' , (12) 

where the uncertainty contains both the uncertainty in (r^) and the systematic uncertainty introduced by 
evaluating the integral only up to finite values (we allow a large range from 1 GeV to 16 GeV for the 
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Fig. 1. The scattering phase in the vector channel (left) for the Breit-Wigner model (dashed line) and the Inverse Amplitude 
Method (solid line). We also plot the square form factor modulus (right). To be able to plot the spacelike and timelike data 
together, the first is plotted against the unphysical variable — \J —fp' with the (negative) spacelike momentum transfer. Here 
we use (r^) as input as described in the text. 

variation of the cut-off, although the integral is basically converged for a cut-off of 2 GeV). It agrees to 
that from vector meson dominance which is about 3.5 GeV~* [2] and it is consistent with the constraint 
[2.3 GeV~'*,5.4 GeV^^] from analyzing the form factor data using analyticity [T^. The advantage of our 
analysis is that it allows in addition for a controlled estimate of the uncertainty. Equivalently, the result in 
term of quartic radius is 

(r-*) = 0.68 ± 0.06 fm'' . (13) 

As mentioned above we will investigate the quark mass dependence of the pion form factor based on the 
assumption that gpjrTr is independent of the quark mass with the dependence of mp taken from other 
sources. Since both parameters are explicit in the parametrization given above, we may study the resulting 
quark mass dependence of Cy, once that of (r^) is fixed. 



3. Chiral perturbation theory 

3.1. General considerations 

In order to determine the quark mass dependence of the square radius, which is the input needed for the 
formalism described above, we will use the results of xPT. Clearly, the curvature Cy as well as its quark 
mass dependence, could also be determined in xPT directly. Depending on the fit and systematics chosen in 
Ref. [16], which is two-flavor 0{p^) xPT calculation, its value could vary between 2 — 6 GeV^"', although the 
authors quote a value around 4 GeV"**, in agreement with a previous estimate [2] (By fitting to form factor 
data, they obtain 3.85 GeV""*). A fit in three-flavor xPT leads to a range 4.49 ± 0.28 GeV""* [13]. 

Adopting = 4± 2 GeV""' as the NNLO xPT resuh, we obtain {r^)/(r'^y = 4± 2. This value is copied 
into Table [II 



3.2. Matching the Omnes representation 



We start by giving the chiral expansion of the vector form factor ^26| valid to NLO in xPT, 



1 , -, , t 1, 



m = 1 + 6^(^ - 4™')^"W + ^^(^6 - 3) . (14) 



(To this order we are free to change M, F to the physical TOtt, fn, since the difference is of NNLO). In this 
expression, 



Jit) 



1 



167r2 



crlog 



cr- 1 
<7 + l 



(15) 



with cr = \/l — Am^/t. A common strategy is to fix the /g constant from the square charge radius [26] 
1 



167r2/- 



■(^6-1) 



(16) 



which is correct up to 0{m'^) in xPT. Higher orders in the chiral expansion cannot bring in powers of t 
since, by definition, the charge squared radius is proportional to the coefficient of the term linear in t in the 
form factor. However, they can bring additional constants to the right hand side (each of a natural order 
suppressed by additional factors of l/(47r/7r)^), and, more important for our purposes, a polynomial of m^. 
To make sure we are not eschewing a critical TO,r dependence, we will compare the right-hand-side of eq. ()16p 
with the NNLO correction in chiral perturbation theory [TB]. The NLO result eq. that depends only 

logarithmically on the pion mass (see eq. ([20]) below), is then extended to 



(0 = 



1 



167r2/2 



1 + i^e - 1) + 



8^2/2 



16^2/2 



167r^ 



13 
192 



181 

Is" 



with 



t2 r 



^6 = ^6 + 6 



f2 



log 



19 
12 



(17) 



(18) 



where ry-^ is a counterterm to be determined empirically, and we will use the simple VMD estimate from 
the same work, at the p scale. 



a("^p) 



-0.25 X 10" 
1.44 



With this estimate, those authors find 

Iq = h 

(the scale-dependence of this number cancels in Eq. (fT8|) . The estimate is taken with constants corresponding 
to set I that we copy in Tabled]). 

Here we have to recall the pion mass dependence of the Vs. The k, as coefficients of the expansion in powers 
of of the Lagrangian density, are by definition pion-mass independent, and so are their renormalized 
counterterms . However, the barred quantities are related to them by absorbing a chiral logarithm 



It 



327r2 



h + log 



(19) 



that makes the Z's scale-independent, but in exchange, pion-mass dependent. This dependence needs to be 
kept track of in the calculation. It becomes crucial in the chiral limit when the pion radius diverges due to 
the virtual pion cloud becoming long-ranged as the pion mass vanishes. This effect appears through /g. 



Therefore, we denote by 



,phys 



the value that the low energy constants take by fitting to physical-world 



data. From here on, when varying the quark (or pion) mass, one needs to change the constant according to 



h = I 



phys 



log 



(20) 



With this we have all the input ready to use Eq. ([7]) also to establish the quark mass dependence of the 
curvature using the Breit-Wigner representation of the phase shifts and the given assumptions on the quark 
mass dependence of both the p mass and coupling. 

However, before we proceed we introduce a method that allows one to estimate the quark mass dependence 
of the p properties directly from the xPT amplitudes evaluated up to a given order, namely the unitarized 
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Fig. 2. Variation of the elastic tttt phase 6\\ with the pion mass. Left: Breit— Wigner model. Note that for m.,r = 'irn^^^ the p 
(held at constant mass) has already crossed below the tttt threshold and is a bound state. Right; Inverse Amplitude Method. 
The resonance stays above the ■kit threshold, its mass having a slight dependence on m^, until rather high pion masses. 



chiral perturbation theory or the inverse ampHtude method (lAM) . The representation we are going to use 
is consistent with NLO chiral perturbation theory at low momentum, and satisfies exact elastic unitarity, 
fitting the pion scattering data up to 1.2 GeV well. The formalism will be introduced in the next subsection. 



3.3. P-wave tttt scattering 



To derive the expression for the lAM one starts with the on-shell tttt scattering amplitude in NLO xPT 
that, for / = 1, is 



Ai[s, t, u) — A{t, s, u) — A{u, t, s) 



(21) 



with 



A{s, t, u) 



1 



+ TTi^ [3 J(s) {s^ ~ ml) + J(t) -u)- 2mlt + Amlu - 2ml) 

1 



p2 Qpi 

+ J{u) {u{u -t)- 2mlu + Amlt ~ 2ml)] + 
5 



967r2/4 



s - 2m^ 



+ M2- 



(s^ + {t- uf) - imlh - 12m2s + 15m^ 



(22) 



The first term can be identified as the leading order, low-energy theorem |27j . but we express it in terms of 
the physical m^, instead of its leading order value M used in the original expression At the meanwhile, 
we keep the independent pion decay constant F. The quantities F and M are related to the physical 
ones via 



F^fAl- 



16^2/2 



= ml\\ 



327r2/2 

The latter expression introduces Z3 into the last line of Eq. (P^ . 

The projection to the spatial p-wave has the usual factor of 1/2 to avoid double-counting quantum states 
by counting all angular configurations with exchanged identical particles 



aii(s) = 



1 1 
32^2 



dcos 0(cos 9)Ai (s, t{s, cos 9),u{s, cos 9)) 



(23) 



One can organize the chiral expansion as 



(24) 



Table 2 

Values of the low energy constants of the NLO SU(2) Chiral Lagrangian. We employ the last row in the calculation. For 
comparison we give several well— known sets. The error refers to the last significant figure. These are determinations based 
on data alone. Several phenomenological and theoretical predictions based on semi-analytical approaches (large Nc, Dyson- 
Schwinger, resonance saturation, etc.) can be found in the literature 31 . 



LEG 


h 


h 


h 


u 


k 


Gasser-Leutwyler I26| 


-2±4 


6± 1.3 


2.9 ± 2.4 


4.3 ±0.9 


16 ± 1 


Dobado et al. 


-0.6 ± 0.9 


6.3 ± 0.5 


2.9 ± 2.4 


4.3 ±0.9 


16 ± 1 


Bijnens et al. set I |16| 


-1.7 


6.1 


2.4 


4.4 ± 0.3 


16 ± 1 


Bijnens et al. set II |16l 


-1.5 


4.5 


2.9 


4.3 




This work 


0.1 ± 1.5 


6± 1.3 


2.9(fix) 


4.3 ±0.9 


16.6 ±0.4 



but the series truncated at any order only satisfies elastic unitarity perturbatively. This is solved, with the 
first two expansion terms, by the Inverse Amplitude Method [28] that reads (suppressing the spin and isospin 
subindices) 



,IAM 



is) 



'■LO 



^hO 



is) 



,NLO 



is) 



(25) 



A Taylor expansion of this amplitude returns NLO xPT as usual for a Pade approximant. However elastic 
unitarity is now exact, and the possibility of a zero of the denominator allows for resonances to appear. 
The associated phase shift 



may be directly employed for the time-like form factor through the Omnes representation. A similar proce- 
dure was taken to calculate the scalar and vector form factors of the pion [29 30 . 

The pion mass dependence of the p meson properties were studied in Rcf. [5j and it was found that gpT^n 
depends only very mildly on the quark mass. In the next section we will investigate the consequences of this 
finding on the pion vector form factor. 

The low energy constants necessary to complete the calculation are fit to the phase shifts data and given 
in Tabled! where they are compared to well-known determinations. Note that with the phase shift data one 
can only determine the difference I2 — h which is about 6 as a result of the fitting [35]. Using Eq. ([7]), the 
curvature can then be obtained as 

= 4.00 ± 0.50 GeV^'^, (26) 

or equivalently, 

(r--*) = 0.73 ± 0.09 fm^ (27) 
The quantity depending solely on the phase shift is 

= 2.13 ±0.42 GeV"''. (28) 
These values are to be considered as our results at the physical pion mass. 

4. Pion-mass dependence needed for lattice extrapolation 

4.1. Mass dependence of the phase shift 

In this section we study the pion mass dependence of the pion vector form factor based on both the 
Breit-Wigner model as well as the amplitudes from the lAM. In the left panel of Fig. [2]we plot the variation 
in the isospin-1 p-wave elastic tttt phase shift Sn with the pion mass in the Breit-Wigner model, where 
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Fig. 3. Dependence of the p mass on the pion mass. Here the p mass is defined as the value of y/s, where the tttt p-wave phase 
shifts cross 90 degrees. Shown are the result from the JAM (solid line) and a fit using Eq. 1291 1 (dashed line) to the lattice data, 
shown as solid dots The lowest point is the physical p mass. 



the physical pion mass is denoted by m^'^^ . For smaU increases in the pion mass, with rup being held fixed 
for illustration, we see how the resonance becomes narrower as the pion threshold approaches. Finally, for 
2r7i7r > rTip, the p becomes bomrd and the phase shift starts at 180 degrees in agreement with Levinson's 
theorem with one bound state. 

Next we consider the lAM. Here what is held constant is the renormalized constants in the chiral 
Lagrangian (/[(/i)) since, as discussed, they are by definition independent of the pion mass. The scale- 
independent V s run logarithmically with the pion mass. This dependence and the explicit pion masses in 
the chiral series bring about a small mTr-dependence of the p mass that puts it just above threshold for 
TO,r — ira^"^^. We plot in the right panel of Fig. [2] the resulting phase as a function of the tttt invariant mass 
for different values of the pion mass. 

The prediction of the pion mass dependence of the p mass resulting from the lAM is shown as the solid 
curve in Fig. [31 The parameters used are li = —0.08, 12 = 5.78, /a = 2.9 and ^4 = 4.3. In this figure we also 
show the results of a recent lattice study |4|. For our comparison we choose this one, for it is the simulation 
where the lowest pion masses are used. To allow for a comparison with recent lattice data, here the p-mass 
is defined as the value of -^s, where the tttt phase shift is 90°. The resulting numbers differ somewhat from 
those corresponding to the real part of the pole position in the second Riemann sheet — the latter definition 
of the mass was used in Ref. [5^. For comparison, we also show the very recent lattice data [1] as the filled 
circles with error bars. The agreement of the lAM with the lattice data is rather satisfying. For later use, 
the lattice data are fitted with an expression derived from an extended version of xPT [36j 

mp = ml + Ciml+C2rnl + Cy,ml\og{^^ . (29) 

The parameter is not included in the fit. It is fixed by the condition that nip = 0.77 GeV at the physical 
pion mass. We find for the p mass in the chiral limit = 0.77 ± 0.1 GeV. The resulting parameters are 

Cl = -0.53 ±0.44 GcV"\ c2=2±lGcV"^ cg = -1 ± 3 GeV"^ (30) 

Using the central values, we get the dashed curve as shown in Fig. [31 Note, the uncertainties in the parameters 
show some correlation, however, since we are here mainly interested in a parametrization of the lattice data, 
we may ignore this observation. Since the pion mass grows faster with the quark mass, eventually the p 
becomes bound (just as the J/tp is under the DD threshold), but this happens for yet larger pion masses. 
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Fig. 4. Pion mass dependence of the mean square radius. We show results based on NLO and NNLO Chiral Perturbation 
theory. The Vy-^ parameter is fixed at its VMD value (its pion-mass dependence contributing at NNNLO is neglected). Likewise 
we plot the mass dependence resulting from a once-subtracted Omnes representation. The data is normalized to the radius for 
physical pion mass. 

4.2. Extrapolation in NLO and NNLO chiral perturbation theory 

Space-like form factors are in principle accessible on a lattice. Since these studies usually employ heavier- 
than-real quarks, the pion mass obtained is also higher than the physical pion mass, and an extrapolation 
is necessary. Another extrapolation to low momentum (due to the finite volume enclosed by the lattice) is 
necessary if the mean square and quartic radii are to be extracted. The mean square radius has indeed been 
studied before [15|34j and extrapolation to physical pion masses taken from chiral perturbation theory. It 
would be interesting to have lattice data at several quark masses to test it. 

Momentum extrapolations to = are, in view of the mean quartic radius, non-linear. In the extraction 
of the mean square radius, the authors of [34] quote a 10% systematic error in the lattice extraction due to 
m^/(l GeV^) xPT errors, and 20% due to g^i„/(l GeV^) momentum extrapolation errors. 

The momentum extrapolation however seems to be avoidable with twisted boundary conditions for the 
fermion fields [TS] , and indeed those authors find 

(r-2)330 Mcv = 0.35(3) fm^, {r^)^^^ mcV = 0.42(3) fm^ . 

where the value at the physical pion mass is obtained with the help of the NLO SU{2) chiral Lagrangian. 

4.3. Chiral extrapolation assisted by the Omnes representation 

We have achieved a representation of the form factor based on the Omnes representation, matched to 
low energy xPT. Since we have relatively good theoretical control of the entire construction, we can now 
extrapolate to unphysical quark (pion) masses. 

The parametrization in Eq. ([7]) requires two pieces of input: the pion scattering p-wave phase shift and the 
mean square radius. For the former we may either use the Breit-Wigner model - together with additional 
assumptions on the p properties - or the lAM, where the quark mass dependence is predicted from NLO 
xPT - higher order pion mass dependencies as they arise from NNLO xPT are not yet studied in this 
framework. 

The square radius has a NLO pion mass dependence caused by the chiral logarithm in Iq. This is a 
major effect for pion masses smaller than physical, towards the chiral limit, but for pion masses higher than 
physical (say the 330 MeV where the lattice data is taken), the term from the NNLO Lagrangian density 
might come to dominate, so we employ this too. Finally, we have an order-of-magnitude countercheck at our 
disposal. By employing a once-subtracted instead of a twice-subtracted Omnes representation, we obtain a 
closed form for the mean square radius in terms of the phase shift 
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Fig. 5. Dependence on the pion mass of the mean quartic radius of the pion (left panel) and Cy (right panel). We show results 
based on the Breit-Wigner model and the Inverse Amplitude Method. The bands correspond to the uncertainties from the 
parameters used (the Z^'s for the JAM and the Ci (c.f. Eq. l|29|l ^ for the Breit-Wigner) as well as from the variation of the 
cut— off. 



' f d.Mf) . (31) 



TT 

4mJ 



All three methods are plotted in Fig. 21 

The results of the pion mass dependence of the quartic radius using the Breit-Wigner model and the lAM 
are plotted in Fig.[Ha) and (c), respectivley, with the dependence of the square radius coming from that 
of le as dictated by Eq. (|^D)) . In the lAM, the pion mass dependence of nip is included intrinsically, while 
in the Breit-Wigner method, it is input from fitting to the recent lattice data 4 as described at the end of 
Subsection 14.11 The bands include the uncertainty from varying the parameters within one sigma and that 
from varying the integration cut-off from 1 to 16 GeV. The uncertainty from the cut-off is the dominant 
one. For comparison, we also plot the result of the Breit-Wigner model with fixed p mass as the dashed 
curve in Fig. [5fc). The dependence is smooth up to the point when the rho becomes stable. Here the curve 
ends. Imposing the mjr-dependent p mass as that given by the lAM, the result for the quartic radius in the 
Breit-Wigner model is shown as the solid curve in Fig. [5fc) . 

We can dispose altogether from the explicit pion mass dependence in (r^) by studying the quantity 
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This constant Cy can of course be also studied on a lattice by itself, although its physical interpretation 
is not transparent. But its mass dependence comes from the phase shift alone (c.f. Eq. ([5])), and is not 
compounded with the mass dependence of the square radius. It is therefore this quantity that allows most 
directly access to the pion mass dependence of the p properties. Our results for this quantity are shown in 
Fig. EJb) and (d) using the lAM and the Breit-Wigner model, respectively. 



5. Summary 



Using the Omnes representation for the pion vector form factor, in this paper we improved the existing 
value for the corresponding curvature using as input only the well known tttt phase shifts in the p-wave as 
well as the pion radius. We find (r'*) — 0.73 ± 0.09 fm^ or equivalently cy = 4.00 ± 0.50 GeV"'' which are 
consistent with the results from NNLO xPT (16113] and recent analysis using analyticity [19]. 

In addition we studied the pion mass dependence of the curvature. A modification of the curvature, called 
Cy in the paper, can be represented solely by the tttt p-wave phase shift. We argued that this quantity 
allows for a clean and model-independent alternative access to the pion mass dependence of the p properties 
and would therefore provide a consistency check of the methods to extract physical parameters from lattice 
simulations. A lattice QCD study of the pion curvature would therefore be of high theoretical interest. We 
also argued that the pion square radius is not well suited for this kind of investigation, since additional, not 
so well controlled, theoretical input would be needed in the analysis. Quantities similar to Cy exist also for 
other form factors, and a study of them from both theoretical and lattice sides would be interesting. 



We would like to thank Stephan Diirr and Jose Pelaez for useful discussions. This work is supported in 
part by grants FPA 2004-02602, 2005-02327, FPA2007-29115-E (Spain), and by the Helmholtz Association 
through funds provided to the virtual institute "Spin and strong QCD" (VH-VI-231). FJLE thanks the 
members of the IKP ( Theorie ) at Forschungszentrum Jiilich for their hospitality during the preparation of 
this work and the Fundacion Flores Valles for economical support. 



References 

[I] M. E. Rose, Phys. Rev. 73 (1948) 279. 

[2] J. Gasser and U.-G. MeiBner, Nucl. Phys. B 357 (1991) 90. 

[3] S. Aoki et al. [CP-PACS Collaboration], Phys. Rev. D 60 (1999) 114508 arXiv;hep-lat/9902018'; S. Diirr et al, Science 
322 (2008) 1224; C. Gattringer, C. Hagen, C. B. Lang, M. Limmer, D. Mohler and A. Schafer, arXiv:081 2.1681l [hep-lat]: 
P. Dimopoulos, C. McNeile, C. Michael, S. Simula and C. Urbach [ETM Collaboration], arXiv:0810.1220 [hep-lat]. 

[4] M. Gockeler, R. Horsley, Y. Nakamura, D. Pleiter, P. E. L. Rakow, G. Schierholz and J. Zanotti, arXiv:0810.5337 [hep-lat]. 

[5] C. Hanhart, J. R. Pelaez and G. Rios, Phys. Rev. Lett. 100 (2008) 152001 arXiv:0801.2871 [hep-ph]]. 

[6] R. R. Akhmetshin et al, JETP Lett. 84 (2006) 413 [Pisma Zh. Eksp. Teor. Fiz. 84 (2006) 491] ar Xiv:hep-ex/0610016] ; 
R. R. Akhmetshin et al. [CMD-2 Collaboration], Phys. Lett. B 648 (2007) 28 arXiv:hep-ex/061002l]; 

[7] A. Aloisio et al. [KLOE Collaboration], Phys. Lett. B 606 (2005) 12 arXiv:hep-ex/0407048 . 

[8] M. N. Achasov et al, J. Exp. Theor. Phys. 101 (2005) 1053 [Zh. Eksp. Teor. Fiz. 101 (2005) 1201] larXiv:hep-ex/0506076] . 

[9] E. P. Solodov [BABAR collaboration], in Proc. of the e+e~ Physics at Intermediate Energies Conference ed. Diego 
Bettoni, In the Proceedings of e-f- e- Physics at Intermediate Energies, SLAC, Stanford, California, 30 Apr - 2 May 2001, 
pp T03 'arXiv:hep-ex/0107027 . 

[10] S. R. Amendolia et al. [NA7 Collaboration], Nucl. Phys. B 277 (1986) 168. 

[II] J. Volmer et al. [The Jefferson Lab F(pi) Collaboration], Phys. Rev. Lett. 86 (2001) 1713 l arXiv:nucl-ex/0010009] . 
[12] G. M. Huber, J. Phys. Conf. Ser. 69 (2007) 012015. 

[13] J. Bijnens and P. Talavera, JHEP 0203 (2002) 046 arXiv:hep-ph/0203049j. 
[14] W. M. Yao et al. [Particle Data Group], J. Phys. G 33 (2006) 1. 

[15] P. A. Boyle et al, JHEP 0807, 112 (2008) arXiv:0804.3971 [hep-lat]]^ 

[16] J. Bijnens, G. Colangelo and P. Talavera, JHEP 9805 (1998) 014 ar Xiv:hep-ph/9805389] . 

[17] P. Masjuan, S. Peris and J. J. Sanz-Cillero, Phys. Rev. D 78 (2008) 074028 larXiv:0807.4893l [hep-ph] . 



12 



[18] I. Caprini, Eur. Phys. J. C 13 (2000) 471 [arXiv:hep-ph /9907227]. 
[19] B. Ananthanarayan and S. Ramanan. iarXiv:0811.0482 [hep-phj. 
[20] R. Omnes, Nuovo Cim. 8 (1958) 316. 

[21] M. Gimeno-Segovia and F. J. Llanes-Estrada, Eur. Phys. J. C 56 (2008) 557 |arXiv:0805."4145l [hep-th]] . 
[22] H. Leutwyler, arXiv:hcp-ph/0212324 

[23] S. J. Brodsky and G. R. Farrar, Phys. Rev. D 11 (1975) 1 309. 

[24] F. Guerrero, Phys. Rev. D 57 (1998) 4136 ^arXiv:hep-ph/9801305] ; F. Guerrero and A. Pich, Phys. Lett. B 412 (1997) 

382 [arXiv:hep-ph/9707347l . 

[25] J. A. Oiler and L. Roca, Phys. Lett. B 651 (2007) 139 [arXiv:0704!00 39' [hep-ph]]. 
[26] J. Gasser and H. Leutwyler, Annals Phys. 158 (1984) 142. 
[27] S. Weinberg, Physica A 96 (1979) 327. 

[28] A. Dobado, M. J. Herrero and T. N. Truong, Phys. Lett. B 235 (1990) 134. 
[29] T. N. Truong, Phys. Rev. Lett. 61 (1988) 2526. 

[30] F. Guerrero and J. A. Oiler, Nucl. Phys. B 537 (1999) 459 [Erratum-ibid. B 602 (2001) 641] arXiv:hep-ph/9805334l . 

[31] G. Ecker, J. Gasser, A. Pich and E. de Rafael, Nucl. Phys. B 321 (1989) 311; F. J. Llanes-Estrada and P. De A. 
Bicudo, Phys. Rev. D 68 (2 003) 094014 arXiv:hep-ph/0306146; ; T. N. Pham and T. N. Truong, Phys. Rev. D 31 
(1985) 3027; F. J. Ynd urain, |arXiv:hep-ph/0212282| M. R. Pennington and J. Portoles, Phys. Lett. B 344 (1995) 399 
farXiv:hep-ph/9409426' . 

[32] A. Dobado and J. R. Pelaez, Phys. Rev. D 56 (1997) 3057 [ arXiv:hep-ph/9604416l . 

[33] A. Dobado et al., Effective Lagrangians for the Standard Model, Springer- Verlag, Berlin Heidelberg 1997. 

[34] T. B. Bunton, F. J. Jiang and B. C. Tiburzi, Phys. Rev. D 74 (2006) 034514 [Erratum-ibid. D 74 (2006) 099902] 

arXiv:hep-lat/0607001 . 

[35] J. F. Donoghue, arXiv: hep -ph /960735 1 1 

[36] P. C. Bruns and U.-G. Meifiner, Eur. Phys. J. C 40 (2005) 97 [arXiv:hep-ph/0411223| . 



13 



